Shiba multiplets due to Mn impurities in MgB 2 
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We study the effect of magnetic Mn ions on the two-band superconductor MgB 2 , and compute 
both the total and spin resolved scanning tunneling spectrum in the vicinity of the magnetic impurity. 
We show that when the internal structure of the Mn ion's d-shell is taken into account, multiple 
Shiba states appear in the spectrum. These multiplets were missed by previous calculations based 
on simplistic models, and their presence could alter significantly the overall interpretation of local 
tunneling spectra for a wide range of superconducting hosts and magnetic impurities. 
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The interaction between a single magnetic impurity 
and the superconducting host reveals fundamental prop- 
erties of both the magnetic ion and the host material. 
This interaction was first studied theoretically, within the 
framework of BCS superconductivity. During the late 
sixties Shiba |]J showed that a magnetic impurity pulls 
down from the continuum states a pair of bound states 
inside the superconducting gap. Indirect indication for 
the presence of finite spectral weight inside the gap of 
an impure superconductor could be inferred from global 
probes of the density of states. However, direct evidence 
for the existence of the so-called Shiba states requires 
an accurate measurement of the local density of states 
near the impurity. Such a measurement became avail- 
able only recently by using high vacuum, low temper- 
ature scanning tunneling spectroscopy (STS). Yazdani 
and his coworkers imaged |2( the local density of states 
around Mn and Gd impurities deposited onto Nb single 
crystals. They found clear evidence for localized states 
in the vicinity of the magnetic impurities, in qualitative 
agreement with Shiba's original findings, and also with 
their own model calculation based on a non-selfconsistent 
solution to Bogoliubov-de Gennes equations. Quanti- 
tative discrepancies, however, are also clearly present, 
especially when comparing the width and spatial de- 
pendence of the resonances to theoretical expectations. 
The presence of magnetic impurity induced bound states 
in a superconductor was turned around and used, both 
theoreticallyH 11 

and experimentally 6] , as an inves- 
tigative tool to probe the unusual ground state of the 
cuprate superconductors. 

In all the studies carried out so far, the direct com- 
parison between theory and experiment is hampered by 
the fact that the theoretical approaches, following Shiba's 
original work, use only a simplified (predominantly clas- 
sical) model to describe the magnetic impurity, and as- 
sume a single spin one-half electron channel that couples 
to the magnetic impurity. Furthermore, the coupling is 



assumed to be in the s-wave channel, and spin-orbit cou- 
pling is generally ignored. In reality, magnetic impurities 
usually have a much more complicated internal structure 
|20j : The magnetic moments are usually due to low-lying 
and crystal-field split d- or /-levels with multiple occu- 
pancy. The aim of this paper is to demonstrate that this 
internal structure has an impact on the structure of the 
Shiba states, and only part of the physics is captured by 
the simplified model used before. In particular, multiple 
channels of charge carriers couple to the magnetic impu- 
rity through channel-dependent coupling. The combina- 
tion of these ingredients generally lead to the appearance 
of multiple pairs of Shiba states. We compute the spatial 
and spin structure of the scanning tunneling microscopy 
(STM) spectra around the magnetic impurity and show 
that these states appear as distinct resonances inside the 
superconducting gap, and can be most clearly resolved 
in spin resolved STM spectra. 

In the following we illustrate our results on the spe- 
cific case of Mn-doped MgB 2 , but we wish to empha- 
size that much of our discussions carry over to other 
systems as well, and that our conclusions are rather gen- 
eral. There are several reasons to choose the Mn — MgB 2 
system. First, in order to observe a Shiba state by scan- 
ning tunneling spectroscopy (STS), one needs a relatively 
large gap. MgB 2 is a perfect candidate in this respect 
since it is a conventional superconductor that has an un- 
usually high critical temperature, T c = 39 K J7j. Sec- 
ond, MgB 2 has a hexagonal AlB^-tyve structure and a 
highly anisotropic band structure [a |9| . As we shall see 
below, this leads to a clear separation of the multiple 
Shiba states. Finally, there is increasing experimental 
evidence from spectroscopic 0, 0, , photoemmision 
[l3| and transport Q measurements that favor the pres- 
ence of two different superconducting gaps in MgB 2 . It 
is therefore an interesting question, how the presence of 
these two gaps influences the structure of Shiba states. 
As mentioned above, MgB 2 crystallizes in the hexagonal 
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AlB2-type structure [j| in which the B~ ions constitute 
graphite-like sheets in the form of honeycomb lattices 
separated by hexagonal layers of Mg ions. Band struc- 
ture calculations 9] indicate that Mg is substantially ion- 
ized, and the bands at the Fermi level derive mainly form 
Boron p orbitals. Four of the six p bands cross the Fermi 
energy, and the Fermi surface consists of quasi-2Z? cylin- 
drical sheets, due to B - p XyV orbitals, and a 31? tubular- 
network (mostly originating from B - p z orbitals). It is 
believed that both structures participate in the forma- 
tion of the superconducting state, though the gap is very 
different on the tubular network and on the cylindrical 
sheets. 
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FIG. 1: Fig. a: Local environment of an Mn ion in MgB 2 . 
The Mn ion is located in a hexagonal cage. Fig. b: Level 
structure and crystal field splitting of the Mn 2+ core states. 



We begin the quantitative analysis of the magnetic im- 
purity problem by establishing how the Mn spins couple 
to the conduction electrons. The interaction part of the 
Hamiltonian depends on the specific location and elec- 
tronic structure of the magnetic impurity considered. In 
what follows, we provide a detailed analysis for Mn im- 
purities, which presumably substitute the Mg ions, and 
take an Mn 2+ c onfig uration with a half-filled d-shell and 
a spin S — 5/2 [21|, though similar considerations hold 
for other types and positions of magnetic impurities. As 
shown in Fig. ^ the five-fold degeneracy of the e?-states 
is lifted by the local hexagonal crystal field into three 
multiplcts that we can label by the original angular mo- 
mentum quantum numbers m of the d-states \m). Each 
of these states is occupied by a single electron, and hy- 
bridizes through a hybridization V m with a specific local 
combination of p-states, tp m , that we constructed explic- 
itly in the framework of a simple tight binding descrip- 
tion of the p-bands 18|. The ipmS transform the same 



way under the hexagonal point group as the |m)'s. In- 
corporating the local Hund's rule coupling and Coulomb 
interaction and the d-states and integrating out virtual 
charge fluctuations to them we arrive at the following 
exchange Hamiltonian: 



= 2 E 

m=-2 



J, 



s (^w™) , (1) 



where the exchange couplings satisfy J m = J_ m ~ 
V^/AE, with AE the characteristic energy of charge 
fluctuations, and the er's denote the Pauli matrices. Note 
that Eq. must conserve m by symmetry, and there- 
fore five independent channels couple simultaneously to 
the impurity. This is the ultimate reason for the appear- 
ance of multiple Shiba states. 

To analyze the Shiba states in detail, we described 
the superconducting state of MgB 2 by the following non- 
interacting BCS Hamiltonian: 



Hn = 



a,k,cr 



a . k 



A a ( 



„t 



h.c. 



Here cj cao . creates an electron with momentum k, spin a, 
energy £k,a, and wave function (j)a,k in one of the four 
bands that cross the Fermi energy (a = 1...4). The 
kinetic energies £k,a and the corresponding wave func- 
tions have been constructed using a simple tight binding 
model, which has proved quite accurate in the vicinity of 
the Fermi surface ^^|. We assume further that the su- 
perconducting gaps take two different values in the p XtV 
and p z bands A xy « 7.5 meV and A z ss 2.5 meV. 

As a next step, we express the interaction part 

J 



Eq. 0} in terms of the operators 



and re-express 



the Hamiltonian in terms of Nambu spinors fl9j <E>k,a = 



{ c k,a,Ti c k,a,J.7 C -k,a,J.' C 



!_ k a |} to obtain 



H = Y,®i a ( £ K a T z + A a T x )$ Ka (2) 

k,a 
k,k' 

The r l 's here denote Pauli matrices acting in the pseu- 
dospin (charge) index of the Nambu spinor. In course 
of the derivation we made use of time reversal symmetry 
and doubled the Hilbert space so that the components of 
the Nambu spinors in Eq. [21 must be considered as inde- 
pendent variables. The form factors / m ,a(k) are normal- 
ized such that they satisfy the orthogonality relation on 
the Fermi surface S a : 



d k/m,a(k.)/ m /_ Q (k) — dm.m' £> a 



(3) 



We determined both the couplings J^ 3 and the form fac- 
tors within the tight binding model of MgB 2 In 
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this approximation - apart from an overall prefactor in 
- they uniquely depend on the crystal structure of 
MgB 2 . The couplings J%f in Eq. approximately sat- 
isfy » (J m a Jm) 1 ^ 2 and are equal in channels ±m 
by symmetry. In the rest of the paper, we neglect quan- 
tum fluctuations of the Mn spin and treat S as a classical 
variable. In this limit we can solve the impurity problem 
exactly. To determine the STM spectrum we first com- 
pute the Nambu Green's function in the presence of the 
impurity: 

G(a,/3,k,k» - 4,k'5 Q ,^(°)(k, W ) (4) 
+ G£°) (k, W ) K WM f m ,^)Gf ( k '> w ) ' 



where G« (k, a>) = (w — e a ^r z — A a r x )~ L denotes the 
non-interacting momentum space Nambu Green's func- 
tion in band a. Eq. Jjyl guarantees that the quantum 
number m is conserved and thus the T-matrix T^ m ' can 
be computed independently for each channel m. The 
computation simplifies further because depends on 
k only through the energy £ Q ,k. As a result, we can ex- 
press T m ^ simply as 



T TO (w) = J m S a/2 [1 - FHJ m S a/2] 



(5) 



where we introduced a matrix notation in the band in- 



dices l a!3 — > T T"" 3 



The matrix F a ^(ui) 



5 a> j3F a (u}) — > F in Eq. OOsimply denotes the local Nambu 
propagator: 



(0) 



F Q (w) = / de g a 



ET 



A Q T* 



with g a the density of states at the Fermi energy in 
band a. For \ui\ < A a the function F a (uj) has only real 
parts and simplifies to F a (ui) = —Trg a (u + A a r x )/(A 2 t — 
w 2 ) -1 / 2 , while for \u\ > A a it is purely imaginary. 

Impurity bound states and resonances can be iden- 
tified from the pole structure of the T-matrix: True 
bound states correspond to zeros of the determinants 
det{T^ 1 (ti>)} on the real axis, and must satisfy \u>\ < A a 
for all bands. Zeros in the vicinity of the real axis, on 
the other hand, correspond to resonances. We found that 
each of the five channels generates a pair of bound states, 
two of which are doubly degenerate by symmetry. 

Let us first discuss the case of two bands only. There 
the calculations further simplify and the energies E of 
the bound states satisfy the following implicit equations: 



{l-g 1 n \a l {±E)){l-gf l a 2 {±E)) 
(g%) 2 ai(±E)a 2 (±E) , 



(7) 



with g^ = TrS^/Qa gpJ^f/2 the dimensionless couplings 
in channel m, and 



a a (u>) = 



A a 



A a 



1/2 



In the absence of inter-band coupling, g 12 = Eq. Q 
would have two independent solutions for each m, and 
the presence of the second superconducting gap would 
double the number of Shiba states. For small g 12, s the 
bound states with energies A z < \E\ < A xy evolve into 
resonances due to mixing with the continuum. We have 
in fact shown that in reality g 12 is large. Consequently, 
the naive expectation mentioned above - based solely on 
the straightforward extension of Shiba's result for two 
bands - is not correct. Using the approximate relation 
for {g 12 f 



9m 9 m> we obtain that 



g 1 m l a 1 (±E)+g%a 2 (±E) = l, 



(8) 



which has only one pair of solutions for each m. Exchange 
coupling between the two bands thus removes half of the 
resonances, and the number of Shiba states is not in- 
creased due to the presence of the second gap in this case. 
We thus obtain altogether five pairs of Shiba states cor- 
responding to the five channels, but two of them are two- 
fold degenerate because of the symmetry = J_ m . 

We are most interested in the total and spin resolved 
tunneling density of states in the vicinity of the Mn sites, 
which are directly measured by conventional and spin po- 
larized STM techniques [22]| . respectively, and are related 
to the local Nambu Green's function G(R,p Q ,cj) at site 
R and orbital p a (a — x, y, z) as 



c ,a(w) = ~ j-ImTr j G(R, p a , u) 1 + Tz 

Z7T Z 



g s , a {uj) = -— ImTr |g(R,p q , uj) crn 
Zn (_ 



1 + T Z 



(9) 



(10) 



with n a unit vector pointing in the direction of the Mn 
spin. These can be computed easily from Q in terms of 
the electronic wave functions 4> a ,\n fis| . 

To obtain a quantitative estimate for the STM spec- 
tra we performed a lengthy, but straigthforward tight- 
binding calculation to determine numerically the form 
factors, the exchange couplings and the electronic wave 
functions above in various geometries. In particular, we 
generalized the discussions above for the case of a semi- 
infinite half plane and computed the STM spectra as a 
function of the Mn positions 01- I n these calculations 
we assumed that the value of the order parameter does 
not depend on the distance from the surface. 

In Fig. [2] we show the density of states for the case 
of a half-infinite sample for various values of the dimen- 
sionless exchange coupling, g = | J2 m a feC) with g a 
the density of states in band a per one spin direction 
at the Fermi energy. We assumed that MgB 2 is cleaved 
along a Boron sheet and that the Mn ion sits right be- 
low this Boron layer. In the STM experiment one pre- 
sumably tunnels into the p z orbitals oriented perpendic- 
ularly to the Boron plane. Therefore we plotted the 
DOS that corresponds to a local p z orbital at position 
A. For generic values of the exchange coupling usually 
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FIG. 2: Normalized density of states for site A as function 
of frequency for different values of g. 



two well-separated pairs of resonances can be observed, 
corresponding to the m = ±1 and m — channels. The 
exchange couplings in channels m — ±2 are much smaller 
than those in channels m = ±1 and m = 0, therefore the 
bound state appears at the superconducting coherence 
peak. The wave functions of the Shiba states and thus the 
amplitudes of the corresponding resonances in the spec- 
trum depend a lot on the tunneling position |18|. More- 
over, as we show in Fig. EI the Shiba states are strongly 
spin-polarized. Therefore, even if a peak happens to be 
close to the BCS coherence peaks, a spin-polarized STM 
can distinguish the Shiba state from the coherence peak. 
Thus, spin-polarized STM would provide a perfect tool 
to identify these multiple Shiba states. 

The presence of spin-orbit coupling should not influ- 
ence these results: Since the ±m (i-levels form time re- 
versed pairs, the symmetry J m — J_ m should hold even 
if we take spin-orbit coupling into account. We conjec- 
ture therefore that the Shiba states remain degenerate 
even in the presence of spin-orbit coupling. 

Finally, let us shortly discuss how these results general- 
ize for other compounds. Possible crystal field structures 
in other cases and the corresponding exchange Hamilto- 
nians have been discussed for various point symmetries in 
the seminal paper of Nozieres and Blandin [2fJ . Though 
the form of the exchange Hamiltonian depends a lot on 
the specific material, in many cases, similar to Mn-doped 
MgB 2 , several channels of conduction electrons couple to 
the local impurity degrees of freedom, and result in mul- 
tiple Shiba states. Thus the appearance of multiple Shiba 
states seems to be a rather general phenomenon. 

In conclusion, we have shown that a magnetic impurity 
generally induces multiple Shiba states in the electronic 
structure of MgB 2 . In particular for Mn we found five 
pairs of Shiba states in the gap, two of which were two- 
fold degenerate. We computed both conventional and 
spin-resolved STM spectra near the impurity site and 
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FIG. 3: Spin polarization and spin resolved density of states 
(inset) for site A and g — 0.491. 



showed that these states can be clearly resolved by both 
methods. Similar multiple Shiba states should appear 
in other superconductors due to the internal structure of 
the magnetic impurity. 
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